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ABSTRACT 

Free-flight  tests  were  conducted  in  the  Defence  Research  Establishment 
Valcartier  (DREV)  aeroballistic  range  on  a  full-scale  MPB-HD  (BDU-5003/B  MOD  1) 
bomb  at  subsonic  velocities.  All  the  main  aerodynamic  coefficients  and  dynamic  stability 
derivatives  were  very  well  determined  using  the  six-degree-of-freedom  single-  and 
multiple-fit  data  reduction  techniques.  The  free-flight  data  shows  that  this  bomb  is 
dynamically  unstable  at  low  angles  of  attack.  The  second  order  pitch  damping  coefficient, 
the  yaw  axial  force  term  as  well  as  side  moments  were  also  reduced.  Wind  tunnel.  Open 
Jet  Facility  experimental  results  and  full-scale  aircraft  free-flight  trials  were  compared 
with  the  aeroballistic  ones. 


RESUME 

Des  essais  en  vol  libre  ont  ete  effectues  dans  le  corridor  aerobalistique  du  Centre 
de  recherches  pour  la  defense  Valcartier  (CRDV)  avec  la  bombe  pleine  echelle  MPB-HD 
(BDU-5003/B  MOD  1)  aux  vitesses  subsoniques.  Tous  les  coefficients  aerodynamiques 
principaux  et  les  derives  de  stabilite  dynamique  ont  et 6  tres  bien  determines  avec  les 
methodologies  de  reduction  de  six  degres  de  liberte  par  les  options  de  reduction  simple  et 
multiple.  Les  donnees  en  vol  libre  demontrent  que  cette  bombe  est  dynamiquement 
instable  a  faible  angle  d'incidence.  Le  deuxieme  terme  du  coefficient  d'amortissement  de 
tangage,  le  terme  du  deuxieme  ordre  pour  la  force  axiale  ainsi  que  les  moments  de  cote 
ont  aussi  ete  reduits.  Les  resultats  experimentaux  de  la  soufflerie  a  rafale  et  a  jet  libre 
ainsi  que  des  essais  d’adronefs  sont  compares  a  ceux  obtenus  dans  le  corridor 
aerobalistique. 
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EXECUTIVE  SUMMARY 


The  CF  have  developed  a  Store  Separation  Model  (SSM)  to  predict  the  separation 
of  stores  from  the  CF-18  aircraft  given  a  configuration  and  initial  conditions.  This  model 
was  developed  in  order  to  reduce  the  risks  of  flight  test  incidents,  and  to  reduce  store 
separation  work  by  directing  efforts  to  critical  areas.  SSM  has  been  used  extensively  by 
Canadair  on  behalf  of  DND  to  support  various  CF-18  stores  clearance  projects  in  the 
past.  The  current  flight  matching  technique  uses  a  trial  and  error  approach,  which  is  very 
time-consuming  and  costly.  It  was  shown  recently  that  the  implementation  of  the 
Maximum  Likelihood  Method  (MLM)  in  the  SSM  could  resolve  its  inherent  deficiencies. 
The  MLM  has  the  capability  of  extracting  aerodynamic  coefficients  and  interference 
parameters,  simultaneously  from  measured  store  separation  trajectories.  The  Ballistic 
SSM  (BSSM),  under  development,  would  be  able  to  predict  full-scale  separation  and 
ballistic  flight  test  data  for  the  CF-18  aircraft. 

Even  though  the  MLM  is  a  well-proven  technique  to  extract  interference 
coefficients  and  aerodynamic  coefficients  (static  and  dynamic),  the  store  separation  tests 
usually  do  not  have  enough  angular  and  translational  motion  so  that  it  can  be  utilized  to 
its  maximum  efficiency.  It  is  therefore  required  to  have  a  very  good  free  stream 
aerodynamic  (static  and  dynamic)  coefficient  database  of  stores  dropped  from  the  CF-18 
to  be  able  to  extract  the  interference  coefficients  with  a  high  degree  of  confidence.  If  the 
free  stream  aerodynamics  of  the  store  are  in  error,  the  MLM  will  over-  or  underestimate 
the  interference  coefficients  to  fit  the  overall  observed  motion.  This  reliable  free  stream 
aerodynamic  database  will  also  be  used  with  the  BSSM  and  the  DREV  6DOF  trajectory 
simulation  program  to  predict  accurate  store  impact  at  the  target  and  in  the  CF-18 
Ballistic  Integrator  Algorithm  OFP. 

DREV  has  a  unique  free  -  flight  aeroballistic  range  where  aerodynamic 
coefficients  (static  and  dynamic)  are  reduced  from  measured  trajectories  with  the  MLM 
methodology.  Projectiles  (scaled  or  full  scale)  are  fired  from  a  powder  gun  through  54 
indirect  shadowgraph  stations.  This  aeroballistic  range  has  shown  over  the  years  to  be 
able  to  extract  very  reliable  aerodynamic  coefficients. 

DREV  was  tasked  by  NDHQ  to  fire  a  first  series  of  store  configurations  in  the 
DREV  aeroballistic  range  with  the  goal  of  obtaining  their  free  stream  static  and  dynamic 
aerodynamic  coefficients.  A  second  objective  was  to  evaluate  various  available  tools  to 
predict  the  aerodynamic  coefficients  in  view  of  reducing  costs  of  experimental  tests.  The 
stores  that  were  chosen  for  this  first  phase  were:  LDGP  MK82  CF,  BDU-5002/B  Mod  1 
(Modular  Practice  Bomb  -  Low  Drag),  and  BDU-5003/B  MOD  1  (Modular  Practice 
Bomb  -  High  Drag).  The  Mach  number  range  of  interest  is  between  Mach  0.6  and  1.5. 

This  report  presents  the  aerodynamic  coefficients  and  stability  derivatives  that 
were  deduced  from  free-flight  tests  conducted  in  the  DREV  aeroballistic  range  on  a  full- 
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scale  BDU-5003/B  Mod  1  (Modular  Practice  Bomb  -  High  Drag)  bomb.  All  the  main 
aerodynamic  coefficients  and  dynamic  stability  derivatives  as  well  as  nonlinear  ones 
were  determined  using  the  six-degree-of-freedom  single-  and  multiple-fit  data  reduction 
techniques.  The  measured  angular  motion  of  all  the  shots  fired  in  the  aeroballistic  range 
shows  a  growth  as  the  projectile  flies  downrange.  This  indicates  that  this  bomb  is 
dynamically  unstable  at  low  angles  of  attack.  The  second  order  pitch  damping  coefficient, 
the  yaw  axial  force  term  as  well  as  side  moments  were  also  reduced.  Wind  tunnel.  Open 
Jet  Facility  experimental  results  as  well  as  full  scale  free-trials  were  compared  with  the 
aeroballistic  range  ones. 

The  data  from  these  trials  combined  with  the  ones  from  the  DREV  Open  Jet 
Facility  suggests  that  this  bomb  flies  with  a  limit  cycle  of  approximately  8.0°  and  12.0° 
angle  of  attack. 

This  database  of  aerodynamic  coefficients  generated  by  this  experimental 
program  is  considered  a  success  and  it  can  be  used  with  the  BSSM,  the  SSM  as  well  as 
the  DREV  6DOF  trajectory  program  to  predict  accurate  weapon  performance.  Also,  it 
can  be  exploited  with  confidence  for  the  CF-18  OFP.  The  methodology  should  be 
expanded  to  other  weapon  systems. 
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NOMENCLATURE 


Variable 

Computer  Output 

Description 

A 

Cross  sectional  area  of  projectile  (m^) 

d 

Diameter  of  projectile  (mm) 

c.  g. 

CG 

Center  of  gravity  (m) 

c'p 

Roll  damping  moment  coefficient 

C16 

Roll  moment  coefficient  due  to  fin  cant 

Cjy 

cig 

Induced  roll  moment  coefficient 

Gnp 

Cap 

Magnus  moment  coefficient 

CQy 

Cng 

Induced  yaw  moment  coefficient 

Gnsm 

Cnsm 

Side  moment  coefficient 

CN 

CN 

Normal  force  coefficient 

CN55a 

CNda 

Trim  force  coefficient  component 

GNs8B 

CNdB 

Trim  force  coefficient  component 

gm 

Cm 

Static  pitch  moment  coefficient 

CMq 

Cmq 

Pitch  damping  moment  coefficient 

CM58A 

Cmda 

Trim  moment  coefficient  component 

CM5SB 

CmdB 

Trim  moment  coefficient  component 

Cmy 

Cmg 

Induced  pitching  moment  coefficient 

cxo 

CXO 

Axial  force  coefficient  at  zero  angle  of 

attack 

CYp 

CYp 

Magnus  moment  coefficient 

CYy 

CYg 

Induced  normal  force  coefficient 

CZy 

CZg 

Induced  normal  force  coefficient 
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Ixj  ly 

- 

Axial  and  transverse 

moments  of  inertia  (kg  m2) 

1 

- 

Length  of  projectile  (m) 

1/d 

- 

Length-to-diameter  ratio 

m 

- 

Mass  of  projectile  (kg) 

M 

Mach 

Mach  number 

MPB-HD 

- 

Modular  Practice  Bomb  -  High  Drag 

P 

- 

Spin  rate  (rad/s  or  deg/m) 

Rei 

- 

Reynolds  number  based  on  length  of 
projectile 

U,  V,  w 

- 

Projectile  component  velocities  (m/s) 

V 

- 

Total  projectile  velocity  (m/s) 

X,  Y,  Z 

- 

Projectile  coordinates  (m) 

t 

- 

Time  of  flight  (s) 

a 

a 

Total  angle  of  attack  (deg) 

amax 

AMAX 

Maximum  angle  of  attack  (deg) 

,Lp 

LN,  LP 

Nutation  and  precession  damping  (1/m) 

0,  \\f,  (j) 

- 

Projectile  orientation  (deg) 

5 

- 

Fin  cant  angle  (rad  or  deg) 

Sp 

- 

Total  trim  angle  (rad  or  deg) 

52 

DBSQ 

Mean  squared  yaw  (deg^) 

e 

- 

Sine  of  the  total  angle  of  attack, 

2  ,  2 
_  V  +  W 
sin  a  = - - — 

V2 

P 

- 

Air  density  (kg/m^) 

6D0F 

__ 

Six  degree  of  freedom 
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Examples 
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1.0  INTRODUCTION 

The  CF  has  developed  a  Store  Separation  Model  (SSM)  to  predict  the  separation  of 
stores  from  the  CF-18  aircraft  given  a  configuration  and  initial  conditions  in  order  to 
reduce  the  risks  of  flight  test  incidents,  and  to  reduce  store  separation  work  by  directing 
efforts  to  critical  areas.  SSM  has  been  used  extensively  by  Canadair  on  behalf  of  DND  to 
support  various  CF-18  stores  clearance  projects  in  the  past.  Because  of  inherent  model 
limitations,  it  is  essential  to  implement  the  capability  to  adjust  aerodynamic  coefficients 
from  the  SSM  database  to  match  model  predictions  with  flight  test  data.  The  current  flight 
matching  technique  uses  an  ineffective  trial  and  error  approach,  which  is  very  time- 
consuming  and  costly. 

The  Defense  Research  Establishment  Valcartier  (DREV)  has  successfully 
implemented  a  computerized  system  which  uses  the  Maximum  Likelihood  Method 
(MLM)  to  iteratively  extract  aerodynamic  coefficients  and  interference  parameters, 
simultaneously,  from  the  trajectory  of  test  articles  in  their  aeroballistic  range  and  Open  Jet 
Facility.  The  heart  of  this  system  are  two  computer  programs  known  as  the  Aeroballistic 
Range  Facility  Data  Analysis  System  (ARFDAS,  Ref.  1)  and  Open  Jet  Facility  Data 
Analysis  System  (OJFDAS,  Ref.  2).  OJFDAS,  (Ref.  2),  successfully  showed  that  it  was 
possible  to  extract  store  separation  interference  coefficients  and  free  stream  aerodynamic 
coefficients  (static  and  dynamic),  simultaneously.  Feasibility  work,  which  confirmed  the 
compatibility  of  the  MLM  algorithms  with  the  SSM,  was  carried  out  under  MLM  Phase  1 
efforts  (Ref.  3). 

A  SSM  and  Ballistic  Store  Separation  Model  (BSSM)  compatible  MLM  algorithm, 
known  as  the  Store  Separation  Model  Data  Analysis  System  (SSMDAS),  was  tested, 
under  Phase  1  (Ref.  3)  efforts  and  confirmed  the  ability  of  MLM  techniques  to  correctly 
adjust  aerodynamic  free  stream  and  interference  coefficients  to  match  SSM/BSSM 
predictions  to  full-scale  separation  and  ballistic  flight  test  data  for  the  CF-18  aircraft. 
Implementation  of  such  an  automated  system  will  improve  the  accuracy  and  efficiency  of 
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DND's  SSM  and  BSSM  for  future  store  separation  and  ballistics  work.  Canadair  has 
implemented  the  MLM  m  the  SSM  and  BSSM  (Ref.  4).  The  modified  SSM  and  BSSM 
shall  have  the  capability  of  using  MLM  techniques  to  achieve  a  match  of  SSM/BSSM 
predicted  trajectories  to  actual  observed  separation  and  free-flight  trajectory  data  of  stores 
dropped  from  CF-18  aircraft  during  flight  test. 

Even  though  the  MLM  is  a  well-proven  technique  to  extract  interference 
coefficients  and  aerodynamic  coefficients  (static  and  dynamic),  the  store  separation  tests 
usually  do  not  have  enough  angular  and  translational  motion  so  that  it  can  be  utilized  to  its 
maximum  efficiency.  It  is  therefore  required  to  have  a  very  good  free  stream  aerodynamic 
(static  and  dynamic)  coefficient  database  of  stores  dropped  from  the  CF-18  to  be  able  to 
extract  the  interference  coefficients  with  a  high  degree  of  confidence.  If  the  free  stream 
aerodynamics  of  the  store  are  in  error,  the  MLM  will  overcompensate  for  this,  which 
might  lead  to  errors  in  the  determined  interference  coefficients.  This  reliable  free  stream 
aerodynamic  database  will  also  be  used  with  the  BSSM  to  predict  accurate  store  impact  at 
the  target  and  in  the  CF-18  Ballistic  Integrator  Algorithm  OFP.  An  NRC  report  (Ref.  5) 
also  states  this  requirement  for  a  reliable  aerodynamic  database:  “In  this  component 
approach  to  store  integration,  the  essential  baseline  information  is  the  store  free  stream 
aerodynamics.  The  aircraft  flow  field,  carriage  loads,  and  launch  characteristics  are 
considered  as  interferences  (not  necessarily  small)  to  the  aerodynamic  characteristic  of  the 
store.  Hence,  whether  flight  tests,  ground  tests,  or  computations  are  used,  a  well- 
established  aerodynamic  database  for  the  store  itself  should  be  obtained”. 

DREV  has  a  unique  free-flight  aeroballistic  range  (Ref.  6  and  7)  where  absolute 
aerodynamic  coefficients  (static  and  dynamic)  are  readily  obtainable  from  measured 
trajectories  with  the  MLM  methodology.  Scaled  or  full-scale  projectiles  can  be  fired  from 
a  powder  gun  through  54  indirect  shadowgraph  stations.  This  aeroballistic  range  has 
shown  over  the  years  to  be  able  to  extract  very  reliable  aerodynamic  coefficients. 
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DREV  was  tasked  by  NDHQ  to  fire  a  first  series  of  store  configurations  in  the 
DREV  aeroballistic  range  with  the  goal  of  obtaining  their  free  stream  static  and  dynamic 
aerodynamic  coefficients.  A  second  objective  was  to  evaluate  various  available  tools  to 
predict  the  aerodynamic  coefficients  in  view  of  reducing  the  costs  of  experimental  testing. 
The  stores  that  were  chosen  for  this  first  phase  were:  a  scaled  LDGP  MK-82  CF  (18.6%), 
a  full-scale  BDU-5002/B  Mod  1  (Modular  Practice  Bomb  -  Low  Drag,  MPB-LD)  and  the 
BDU-5003/B  Mod  1  (Modular  Practice  Bomb  -  High  Drag,  MPB  -  HD).  The  Mach 
number  range  of  interest  is  between  Mach  0.6  and  1.5. 

This  report  presents  the  aerodynamic  coefficients  and  stability  derivatives  that 
were  deduced  from  free-flight  trajectories  measured  in  the  DREV  aeroballistic  range  on  a 
full-scale  BDU-5003/B  Mod  1  (Modular  Practice  Bomb  -  High  Drag)  bomb.  All  the  main 
aerodynamic  coefficients  and  dynamic  stability  derivatives  as  well  as  nonlinear  ones  were 
determined  using  the  six-degree-of-freedom  single-  and  multiple-fit  data  reduction 
techniques.  The  measured  angular  motion  of  all  the  shots  fired  in  the  aeroballistic  range 
shows  a  growth  in  the  incidence  as  the  projectile  flies  downrange.  This  indicates  that  this 
bomb  is  dynamically  unstable  at  low  angles  of  attack.  The  second  order  pitch  damping 
coefficient,  the  yaw  axial  force  term  as  well  as  side  moments  were  also  reduced.  Wind 
tunnel.  Open  Jet  Facility  experimental  results  as  well  as  full-scale  free-trials  were 
compared  with  the  aeroballistic  range  ones. 

This  trial  was  performed  at  DREV  in  October  1998  and  the  analysis  in  February 
1999,  under  Work  Unit  3ecl6,  Improvement  to  CF-18  Ballistics  Algorithms. 
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2.0  MODEL  CONFIGURATION 

2.1  MPB-HD  Configuration 

The  full-scale  MPB-HD  was  a  proper  candidate  since  its  caliber  (50.8  mm)  was  a 
suitable  size  to  fire  in  the  DREV  aerobalhstic  range  at  full  scale.  Also,  an  extensive  wind 
tunnel  database  exists  at  DREV  with  this  configuration  and  it  was  also  tested  at  length 
from  aircraft.  This  was  also  one  of  the  configurations  that  were  tested  m  the  DREV  Open 
Jet  Facility  and  used  as  a  basis  to  validate  the  MLM  methodology  m  obtaining  store 
interference  coefficients  and  free  stream  aerodynamic  coefficients,  simultaneously  (Ref. 
2). 


The  in-service  MPB-HD  configuration  was  used.  The  main  dimensions  are 
provided  in  Fig.  1  in  caliber.  The  reference  diameter  is  50.8  mm.  The  only  external 
geometry  difference  from  the  standard  practice  bomb  is  that  the  fms  were  rotated  45°  from 
the  locator  holes,  i.e.  one  pair  of  fins  is  in  line  with  the  locator  holes.  This  was  done  so  as 
to  be  able  to  launch  them  with  the  sabot  that  was  designed  (Ref.  8).  The  fumer  cartridge 
was  also  replaced  by  a  dummy  one  to  keep  the  center  of  gravity  and  the  mass  as  close  as 
possible  to  the  m  service  bomb. 

The  MPB-HD  has  a  1.35  cal  ogive  nose  followed  by  a  4.07  cal  cylindrical  portion 
and  the  fins  are  placed  at  the  end  of  an  2.48  cal  extended  boattail.  The  fms  have  a  2.00  cal 
span  and  are  of  a  clipped  delta  type.  The  fin  leading  edges  are  blunt  with  a  thickness  of 
0.08  cal  which  reduces  to  0.05  cal  at  the  trailing  edge.  The  MPB-HD  has  the  same  fin  and 
body  dimensions  as  the  LD  version.  A  high  drag  0.07  cal  thick  retardation  disk  with  a 
diameter  of  1.76  cal  is  located  at  1.89  cal  from  the  nose.  A  1.7  cal  diameter  high  drag 
conical  tail  is  placed  just  aft  of  the  fms.  The  fin  type  and  main  dimensions  are  the  same  as 
in  the  low  drag  version.  The  fins  have  no  cant  to  produce  spin  rate.  The  center  of  gravity 
of  the  tested  projectiles  was  located  at  3.97  caliber  from  the  nose  with  the  dummy 
cartridge  in  the  projectile.  The  total  length  of  the  projectile  is  8.55  cal. 
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The  nominal  physical  properties  of  the  model  are  given  in  Table  I  and  the  physical 
properties  of  each  test  projectile  are  provided  in  Table  II. 


2.2  Sabot  Design 

Since  the  projectile  had  to  be  launched  from  a  powder  gun  to  conduct  tests  in  the 
DREV  aeroballistic  range,  special  sabots  had  to  be  designed  to  fire  them.  Since  the  model 
configuration  in  this  case  is  fin  stabilized,  a  smooth  bore  gun  was  utilized.  The  standard 
gun  employed  at  DREV  to  fire  fin-stabilized  projectiles  of  these  dimensions  in  the 
aeroballistic  range  is  a  110-mm  smooth  bore  gun. 

Several  aspects  have  to  be  considered  when  designing  sabots  and  models.  They 
are:  projectile  configuration,  total  mass,  sabot  separation  at  the  sabot  trap  located  at  9.2  m 
from  the  muzzle  at  the  aeroballistic  range,  muzzle  velocity  desired,  gun  accelerations,  etc. 
The  last  three  mentioned  have  to  be  consistent  from  round  to  round.  In  these  tests,  the 
highest  muzzle  velocity  achieved  was  approximately  323  m/s  (Mach  0.95)  and  the  lowest, 
270  m/s  (Mach  0.8). 

The  tail  portion  of  the  modular  practice  bombs  is  made  of  a  polycarbonate 
material.  Therefore,  it  is  impossible  to  launch  them  with  a  base  plate  pusher  sabot  since 
the  projectile  would  disintegrate  at  launch.  The  modular  practice  bombs  have  two  locator 
holes  situated  close  to  the  center  of  gravity  (Ref.8)  and  a  sabot  design  that  would  pull  the 
projectile  by  these  holes  presented  an  interesting  option. 

Figure  2  shows  a  schematic  of  the  sabot  with  the  MPB-HD  projectile  inside  it.  The 
sabot  for  the  MPB-LD  is  exactly  the  same.  The  detail  drawings  of  the  sabot  are  provided 
in  Ref.  8.  It  is  a  two  petal  sabot  design  made  of  aluminum.  The  lengths  of  the  saw  cuts  on 
each  side  were  adjusted  to  obtain  adequate  petal  separation  for  the  expected  velocities.  A 
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sabot  base  seal  pad  was  also  used  to  prevent  gas  leakage  past  the  sabot  body.  It  has  two 
pins  at  the  front  of  the  sabot  to  pull  the  projectile  down  the  barrel.  These  pins  were 
designed  to  fit  the  in-service  MPB  locator  holes. 

A  pivot  pin,  which  is  in  line  with  the  saw  cuts,  was  added  to  force  the  sabot 
opening  at  that  point.  A  polycarbonate  ring  with  a  5°  angle  is  positioned  at  the  aft  end  of 
the  sabot.  There  are  two  reasons  for  this.  The  first  one,  is  to  have  a  good  pressure  seal 
between  the  sabot  and  the  gun  tube  so  as  to  be  able  to  have  a  known  shot  start  pressure 
which  helps  in  having  consistent  muzzle  velocities  at  the  same  propellant  charge  mass. 
The  second  reason  is  that,  as  the  sabot  leaves  the  gun  tube,  the  high  radial  pressure  acting 
on  the  rear  ring  relative  to  the  front  part,  causes  the  pivoting  action  at  the  pivot  point  of  the 
sabot  petals. 

A  photograph  of  the  sabot-model  package  as  well  as  all  the  components  is  shown 
in  Fig.  3.  One  should  notice  the  roll  pin  placed  on  one  of  the  model  fins  to  make  it 
possible  to  measure  the  roll  orientation  of  the  projectile  when  fired  in  the  aeroballistic 
range.  The  total  model-sabot  mass  is  approximately  5.1  kg.  The  complete  drawings  of  the 
sabot  as  well  as  the  details  of  the  proof  trials  that  were  conducted  to  verify  sabot-model 
integrity  at  launch  can  be  found  in  Ref.  8. 


3.0  FACILITY  DESCRIPTION  AND  TEST  CONDITIONS 

3.1  DREV  Aeroballistic  Range 

The  DREV  aeroballistic  range  (Refs.  6  and  7)  is  an  insulated  steel-clad  concrete 
structure  used  to  study  the  exterior  ballistics  of  various  free-flight  configurations.  The 
range  complex  consists  of  a  gun  bay,  control  room  and  the  instrumented  range  (Fig.  4a).  A 
massive  blast  wall  is  located  in  front  of  the  building  to  stop  sabot  pieces  and  minimize 
vibrations  transmitted  to  the  range  structure  and  instrumentation.  Projectiles  of  calibers 
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ranging  from  5.56  to  155  mm,  including  tracer  types,  may  be  launched.  Large-caliber 
models  have  been  fired  up  to  Mach  7. 

The  230-meter  instrumented  length  of  the  range  has  a  6.1-m  square  cross  section 
with  a  possibility  of  54  instrumented  sites  along  the  range  (Fig.  4b).  For  these  tests,  41  of 
the  stations  were  utilized.  These  sites  house  fully  instrumented  orthogonal  shadowgraph 
stations  that  yield  photographs  of  the  shadow  of  the  projectile  as  it  flies  down  the  range. 
The  maximum  shadowgraph  window,  an  imaginary  circle  within  which  a  projectile  will 
cast  a  shadow  on  both  reflective  screens,  is  1.6  meters  in  diameter.  There  are  also  four 
Schlieren  stations  (two  operational  for  these  tests)  at  the  beginning  of  the  range  that  yield 
high  quality  flow  photographs.  The  range  is  also  air  conditioned  to  maintain  a  constant 
relative  humidity  of  approximately  45%.  The  nominal  operational  conditions  of  the  range 
are  20°  C  at  standard  atmospheric  conditions.  The  spark  source  and  reference  point 
locations  that  were  used  were  deduced  from  a  standard  survey.  A  dynamic  calibration  was 
conducted  in  the  X,  Y,  X,  0  and  \|/  coordinates. 


3.2  Test  Conditions  and  Particularities 


Eleven  (11)  projectiles  were  fired  in  the  aeroballistic  range  program  with  the  110 
mm  smooth  bore  gun  with  the  HI-LO  adapter  (Ref.  8).  All  the  projectiles  had  roll  pins. 
The  range  conditions  for  each  test  projectile  at  time  of  firing  are  indicated  in  Table  HI. 
The  muzzle  velocities  range  from  a  low  of  270  m/s  (Mach  0.8)  to  a  maximum  of  323  m/s 
(Mach  0.95).  The  mid-range  Mach  numbers  varied  from  0.68  to  0.83  which  yielded 
Reynolds  number,  based  on  the  length  of  projectile,  between  6.6xl06  and  8.0xl06, 
respectively.  The  initial  angles  of  attack  ranged  from  a  low  of  1.5°  to  a  maximum  of  4.5°. 

The  gun  muzzle  was  situated  at  a  downrange  coordinate  of  approximately  6.32  m 
in  the  aeroballistic  range  coordinate  system.  Due  to  the  low  muzzle  velocities,  the 
movable  butt  was  utilized  to  capture  the  projectiles  before  the  end  of  the  range.  It  was 
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placed  after  the  47th  shadowgraph  stations  at  approximately  150.0  m.  This  allowed  a 
possibility  of  41  shadowgraph  stations,  or  115.0  m,  to  measure  data.  This  butt  was  used 
for  all  the  projectiles. 

Due  to  the  length  of  the  models,  the  total  shadow  of  the  projectiles  in  flight  could 
not  be  captured  on  the  films.  A  time  delay  was  deliberately  set  so  that  at  least  the  base  of 
the  model  could  be  seen  with  certainty.  This  allowed  the  reading  of  the  films  at  the  back 
of  the  collar  and  at  the  front  of  the  fins  for  calculating  the  trajectory. 

A  typical  Schlieren  photograph  showing  the  complex  flow  field  and  shock 
structure  of  a  projectile  in  flight  can  be  seen  in  Fig.  5  for  shot  A12  at  Mach  0.94. 

The  numbering  scheme  to  refer  to  the  shots  and  a  particular  configuration  is  as 
follows.  One  letter  followed  by  6  digits,  as  for  example  A981002,  identifies  the  shot 
numbers.  The  letter  corresponds  to  a  particular  configuration.  The  first  four  numbers 
(9810)  indicate  the  date  (year  and  month)  that  the  projectile  was  fired  in  the  range.  The 
last  two  digits  correspond  to  the  shot  number  for  that  particular  configuration.  For  the 
example  given  above,  the  shot  number  corresponds  to  the  second  shot  of  the  Model  A 
configuration  that  was  fired  in  the  range  in  October  1998.  For  convenience,  the  shot 
numbers  are  usually  referred  to  the  letter  and  the  shot  number,  A02. 


4.0  FREE-FLIGHT  DATA  REDUCTION 

Extraction  of  the  aerodynamic  coefficients  and  stability  derivatives  is  the  primary 
goal  in  analyzing  the  trajectories  measured  in  the  DREV  aeroballistic  range.  This  was 
performed  by  means  of  the  Aeroballistic  Range  Facility  Data  Analysis  System  (ARFDAS, 
Ref.  1)  described  in  Fig.  6.  These  programs  incorporate  a  standard  linear  theory  and  a  six- 
degree-of-freedom  (6DOF)  numerical  integration  technique.  The  6DOF  routine 
incorporates  the  Maximum  Likelihood  Method  (MLM)  to  match  the  theoretical  trajectory 
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with  the  experimentally  measured  trajectory.  The  MLM  is  an  iterative  procedure  that 
adjusts  the  aerodynamic  coefficients  to  maximize  a  likelihood  function.  The  application  of 
this  likelihood  function  eliminates  the  inherent  assumption  in  least  square  theory  that  the 
magnitude  of  the  measurement  noise  must  be  consistent  between  parameters  (irrespective 
of  units).  In  general,  the  aerodynamic  coefficients  are  nonlinear  functions  of  angle  of 
attack,  Mach  number  and  roll  angle. 

ARFDAS  represents  a  complete  ballistic  range  data  reduction  system  capable  of 
analyzing  both  symmetric  and  asymmetric  models.  The  essential  steps  of  the  data 
reduction  system  are  to  (I)  assemble  the  dynamic  data  (time,  position,  angles),  model 
measured  physical  properties  and  atmospheric  conditions,  (2)  perform  linear  theory 
analysis,  and  (3)  perform  6DOF  analysis. 

These  three  steps  have  been  integrated  into  data  analysis  system  to  provide  the  test 
scientist  with  a  convenient  and  efficient  means  of  interaction.  At  each  step  in  the  analysis, 
permanent  records  for  each  shot  are  maintained  so  that  subsequent  analyses  with  data 
modification  are  much  faster. 

The  6DOF  data  reduction  system  can  also  simultaneously  fit  multiple  data  sets  (up 
to  five)  to  a  common  set  of  aerodynamics.  Using  this  multiple-fit  approach,  a  more 
complete  range  of  angle  of  attack  and  roll  orientation  combinations  is  available  for 
analysis  than  would  be  available  from  a  single  flight.  This  increases  the  accuracy  of  the 
determined  aerodynamic  coefficients  over  the  entire  range  of  angles  of  attack  and  roll 
orientations. 

The  aerodynamic  data  presented  in  this  report  were  obtained  using  the  fixed-plane 
6DOF  analysis  (MLMFXPL)  with  both  the  single-  and  multiple-fit  data  correlation 
techniques.  The  equations  of  motion  have  been  derived  in  a  fixed-plane  coordinate  system 
with  Coriolis  effects  included.  The  formal  derivation  of  the  fixed-plane  model  is  given  in 
Ref.  9. 
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All  the  results  given  here  were  reduced  after  the  dynamic  calibration  biases  were 
accounted  for  the  X,  Y,  X,  0  and  \j/  coordinates.  The  methodology  of  the  dynamic 
calibration  for  the  DREV  aeroballistic  range  is  explained  m  Ref.  10. 


5.0  FREE-FLIGHT  RESULTS  AND  DISCUSSIONS 


The  aerodynamic  coefficients  and  stability  derivatives  that  were  reduced  from  the 
free-flight  trajectories  measured  in  the  aeroballistic  range  are  presented  in  tabular  form  for 
the  linear  theory  analysis  and  in  both  tabular  and  plotted  format  for  the  6DOF  reductions. 
All  of  the  determined  aerodynamic  coefficients  are  given  at  the  mid  range  measured  Mach 
number. 


5.1  Linear  Theory  Results 

The  linear  theory  parameters  deduced  from  the  decoupled  motion  are  provided  in 
Table  IV.  The  magnitudes  of  the  initial  angles  of  attack  varied  from  a  low  of  1 .6°  to  a  high 
of  roughly  4.5°.  The  amplitude  of  the  initial  nutation  and  precession  arms,  KF  and  KS, 
and  the  mean  squared  yaw  (Dbsq)  provides  an  indication  of  these  angles  of  attack. 

In  all  cases  the  shots  were  dynamically  unstable,  as  observed  by  the  positive 
nutation  and  precession  damping  modes  (LF  and  LS).  This  implies  that  the  angular  motion 
is  increasing  as  the  projectile  flies  downrange.  The  frequencies  (WF  and  WS)  are 
consistent.  It  should  be  noted  that  the  trim  angles  (KT)  are  of  the  order  of  0.5°  in  some 
cases.  This  will  be  further  investigated  in  the  6DOF  analysis. 

The  aerodynamic  coefficients  deduced  from  the  linear  theory  parameters  are 
presented  in  Table  V.  The  methodology  to  obtain  the  aerodynamic  coefficients  is 
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explained  in  Ref.  1.  The  main  aerodynamic  coefficients  (Cxo,  C^a  >  CMa,  Cjp  and 
C1g8)  are  consistent.  The  pitch  damping  term,  C^q ,  is  positive  indicating  a  dynamic 
instability  for  all  the  shots. 

The  standard  deviation  error  in  the  angular  motion  (E-Ang)  from  the  linear  theory 
analysis  (Table  V)  is  high  in  some  cases.  The  dynamic  instability  is  not  well  resolved  with 
the  linear  theory  analysis  and  this  suggests  that  the  linear  theory  analysis  was  not  fitting 
some  parameters  adequately,  probably  due  to  nonlinear  variation  with  angle  of  attack  in 
some  aerodynamic  coefficients.  These  nonlinearities,  if  they  exist,  are  best  modeled  and 
reduced  with  the  6DOF  reduction  technique  of  the  next  section. 


5.2  Six-Degree-of-Freedom  Results 

The  determined  aerodynamic  coefficients,  their  standard  deviation  errors,  and  the 
standard  deviation  errors  between  the  theoretical  and  experimental  trajectories  for  the 
axial,  angular  and  roll  motions  are  given  in  Tables  VI  and  VII  for  the  single-  and  multiple- 
fit  data  reduction  techniques,  respectively.  The  moment  reference  center  for  the  pitch  and 
moment  coefficients  was  at  46.0%  of  the  length  from  the  nose  of  the  projectile  (3.93  cal). 
All  the  results  are  given  at  the  mid-range  Mach  number  for  the  single-fit  data  reductions 
and  at  the  average  mid-range  Mach  numbers  for  the  multiple-fit  data  reductions. 

A  coefficient  that  appears  with  a  value  and  a  (*)  between  parentheses  directly 
below,  indicates  that  this  coefficient  was  held  constant  and  one  that  has  a  (-)  between 
parentheses  indicates  that  this  coefficient  was  solved  for  and  that  the  standard  deviation 
error  for  this  coefficient  was  higher  than  100%,  that  is,  it  does  not  influence  the  fit  and  is 
considered  undetermined.  Those  with  numbers  between  parentheses  represent  the  standard 
deviation  error  for  that  particular  coefficient. 
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The  multiple-fit  groups  were  chosen  by  Mach  numbers  and  three  groups  of 
multiple-fit  data  reductions  were  conducted,  as  given  in  Table  VII.  In  some  instances,  the 
results  showed  more  variation  in  Cxo  and  C^a  than  would  be  expected  at  the  same 

Mach  number.  Therefore,  for  some  multiple-shot  groups,  the  multiple-fit  data  reductions 
were  conducted  as  follows.  The  axial  force  coefficient  was  kept  constant  at  the  average  of 
the  single-fit  results.  The  variation  from  this  value  was  then  uniquely  solved  for  each  shot. 
A  unique  CMa  was  solved  for  each  individual  shot  in  the  cases  where  the  variation  was 

deemed  too  high.  It  is  believed  that  these  variations  originate  from  the  flow  over  the  fins 
due  to  the  high  drag  collar  located  at  the  front  of  the  model.  These  would  be  dependent  on 
the  angle  of  attack  and  roll  orientation  of  the  projectile. 


As  seen  from  Tables  VI  and  VII,  all  of  the  main  aerodynamic  coefficients  (Cxo, 

Cncc’  Cm<x  ar>d  CjyjqQ )  were  very  well  determined,  as  indicated  by  the  low  probable 

errors  of  fits  on  the  coefficients.  It  should  be  noticed  that  the  Cj^qo  values  are  positive. 

The  pitch  damping  quadratic  coefficient  expansion  term  (Cm  )  was  well  determined 

qa2 

in  the  multiple-fit  data  reductions  at  Mach  0.8  and  0.82.  The  single  fit  was  held 

qa2 

constant  at  the  determined  multiple  values  and  Cp^qO  was  solved  for. 


Since  there  was  no  fin  cant  on  the  projectiles,  the  roll  motion  was  very  limited. 
Therefore,  the  roll  damping  moment,  Cjp ,  was  kept  constant  at  the  PRODAS  (Ref.  11) 

estimates  obtained  for  the  MPB-LD  bomb  (Ref.  12)  and  the  roll  producing  moment  due  to 
fin  cant,  Cs§S ,  was  allowed  to  vary  to  take  into  account  any  manufacturing  tolerances  in 
the  fin  angles.  The  aerodynamic  trims  for  the  pitch  moment  and  the  normal  force  were 
solved  for  all  the  shots. 
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The  yaw  axial  force  (Cx  0  )  was  also  well  resolved  in  the  multiple-fit  data 

az 

reduction  and  the  respective  single  fits  were  held  constant  at  the  multiple-fit  determined 
value.  A  pure  side  moment  (Cnsm )  was  resolved  in  two  multiple  shot  groups.  There  is  no 

doubt  that  this  side  moment  originates  from  the  severe  turbulence  caused  by  the  collar  on 
the  fins.  The  projectile  roll  motion  was  almost  nothing  in  most  cases. 

The  standard  deviation  errors  of  the  single  and  multiple  fits  are  of  the  order  of  2.5 
mm  in  the  downrange  coordinate,  1.0  mm  in  the  swerve  motion,  0.6°  in  pitch  and  yaw  and 
of  the  order  of  5.0°  in  roll.  These  errors  of  the  fits  in  pitch  and  yaw  and  in  the  swerve 
direction  are  a  bit  high  when  compared  with  other  test  programs  conducted  in  the  DREV 
aeroballistic  range.  The  aerodynamic  coefficient  expansions  utilized  is  probably  not  as 
well  suited  to  analyze  these  type  of  flows.  The  6DOF  probable  errors  of  fits  are  a  bit 
smaller  than  the  linear  theory  ones  because  of  the  better  mathematical  modeling  of  the 
motion,  such  as  the  inclusion  of  aerodynamic  trims,  angle  of  attack  dependent  terms  and 
variations  with  Mach  number. 


5.3  Comparison  of  6DOF  Single-  and  Multiple-Fit  Results 

A  comparison  of  the  reduced  aerodynamic  coefficients  from  the  6DOF  data 
reductions  techniques  with  the  single-  and  multiple-fit  results  are  given  in  Figures  7  to  13. 
The  single  fit  data  points  (AB  -  SF)  are  shown  as  open  triangles  while  the  multiple-fit  data 
reduction  results  (AB  -  MF)  are  given  as  solid  triangles. 

Appendix  A  presents,  for  every  test  shot,  the  total  angle  of  attack  history  with  the 
observed  angular  motion  and  the  theoretical  determined  one  with  the  reduced  aerodynamic 
coefficients.  The  experimental  data  points  (open  circles)  and  the  calculated  trajectory 
(continuous  line)  from  the  determined  coefficients  are  compared.  This  allows  a 
verification  that  the  reduced  aerodynamic  coefficients  do  fit  the  experimental  trajectory 
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satisfactorily.  For  every  shot,  the  total  angle  of  attack  and  the  angular  motion  plots  in  pitch 
and  yaw  are  given  as  a  function  of  the  downrange  coordinate. 

The  axial  force  coefficient  at  zero  angle  of  attack  (Cxo)  as  a  function  of  Mach 

number  is  shown  in  Fig.  7a.  Cxo  is  of  the  order  of  3.5  with  a  high  scatter  m  the  single-fit 

results.  The  2nd  order  axial  force  coefficient  term,  Cx  -  ,  was  well  determined  and  Cx 

a/ 

as  a  function  of  angle  of  attack  for  the  three  groups  of  multiple  fits  is  shown  in  Fig.  7b. 

An  analysis  was  conducted  for  every  shot  to  investigate  trends  with  angles  of 
attack.  The  data  was  fit  by  sections  (groups  of  15  stations)  for  different  lengths  along  the 
range  as  provided  below: 


1st  station 

Last  station 

1 

15 

7 

21 

13 

27 

19 

33 

25 

39 

The  total  axial  force  coefficient  as  a  function  of  the  mean  angle  of  attack  is  shown  in  Fig. 
7c.  The  variation  of  Cx  as  a  function  of  the  mean  angle  of  attack  is  linear.  The  total  Cx 

determined  for  the  multiple  fits  as  well  as  an  attempt  to  a  fourth  order  expansion  (Cx  4 

=  -10000.0)  is  superimposed  on  the  sectional  data.  The  fourth  order  expansions  term 
appears  to  better  fit  the  data,  but  the  aeroballistic  range  observed  angles  of  attack  were  not 
high  enough  to  confirm  this  hypothesis. 

CNa,  the  normal  coefficient  slope  versus  Mach  number  is  displayed  in  Fig.  8. 
There  is  a  slight  scatter  in  the  single-fit  results  and  C^  is  about  15.0  over  the  whole 
Mach  number  range  tested. 


P515888.PDF  [Page:  27  of  72] 


UNCLASSIFIED 

15 


The  variation  of  the  pitching  moment  coefficient  slope,  CMa ,  with  Mach  number 
is  shown  in  Fig.  9.  There  is  only  a  very  slight  scatter  in  the  single-fit  results  for  CMtx  and 
it  is  roughly  -50.0  between  Mach  0.65  and  0.82. 

The  total  pitch  damping  coefficient  is  defined  as 

CMq  =  CMqO  +  £2 

The  variation  of  the  pitch  damping  moment  coefficient  at  zero  angle  of  attack,  Cj^qO  >  anc* 

the  second  order  expansion  term,  ,  with  Mach  number  are  shown  in  Fig.  10a  and 

qa2 

Fig.  10b,  respectively.  As  explained  previously,  Cp^qO  is  positive  since  the  angular 

motion  increases  with  range  and  the  very  high  negative  values  of  controls  the 

qa2 

amplitude  of  the  limit  cycle.  The  variation  in  the  numbers  are  quite  high  since  no  final 
limit  cycles  were  achieved  in  the  short  range  of  the  tests  and  the  angular  motion  keeps 
increasing. 

Cjyfq  as  a  function  of  angle  of  attack  for  the  three  multiple  data  reductions  are 

shown  in  Fig.  10c.  As  the  angle  of  attack  increases  CMq  crosses  the  zero  barrier  at 

roughly  3.0°.  The  angular  motion  can  best  be  explained  with  the  linear  theory  analysis 
formulation  (Ref.  11).  The  pitch  and  yaw  motions  are  modeled  by  a  damped  sinusoidal 
function.  The  yaw  damping  factors  of  the  nutation  and  precession  arms  for  a  non-rolling 
projectile  with  no  side  moments  and  induced  moment  terms  are  given  by: 
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P  A 
4  m 


-C 


Na 


+  2Cd  + 


kr 


C-Mq 


where 


For  low  drag  projectiles  the  (2Cq  )  term  is  often  neglected  since  it  is  small  compared  with 
.  This  is  not  the  case  for  this  high  drag  configuration. 

The  criterion  for  a  projectile  to  be  dynamically  stable  is  that  the  fast  and  slow  arms 
damping  factors  be  negative.  Dynamic  stability  boundaries  for  the  fast  and  slow  arms 
(which  are  the  same  m  this  case)  for  the  pitch  damping  coefficient  can  be  calculated  from 
the  above  as: 


CMq  (*N  =  0.0  ,  =  0.0)  =  2  [  CNa  -  2Cd]  k? 

Using  Cp>  =  3.5  and  CNot  =  15.0,  yields  a  stability  bound  for  CMq  of  roughly  39.0.  If 

this  stability  bound  is  superimposed  on  Fig.  10c,  it  can  be  seen  that  the  crossover  point  is 
of  the  order  of  3.0°  to  4.0°.  This  agrees  quite  well  with  the  motion  plots  at  the  mid-range 
value. 


The  total  pitch  damping  as  a  function  of  the  mean  angle  of  attack  is  shown  in  Fig. 
lOd.  There  is  no  apparent  pattern  of  CMq  as  a  function  of  the  mean  angle  of  attack  for  the 


different  sections. 
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Ci  ,  the  roll  damping  coefficient,  is  demonstrated  versus  Mach  number  in  Fig.  11. 

ir 

This  term  was  not  solved  for  due  to  the  lack  of  roll  motion.  These  are  the  same  values  as 
the  MPB-LD  (Ref.  12),  which  were  obtained  from  the  predicted  values  from  PROD  AS 
(Ref.  13)  and  it  is  about  -3.9. 

The  trim  moment  and  trim  normal  force  coefficients  were  well  determined  (Table 
VI)  and  they  improved  the  fit  noticeably  when  they  were  included  in  the  fitting  process. 
The  asymmetries  caused  by  modular  design  of  the  MPB-HD  and  the  manufacturing 
tolerances  contribute  mostly  to  these  coefficients.  A  trim  angle  can  be  calculated 
separately  from  the  force  and  the  moment  term.  The  total  moment  trim  is: 

cMg5T  =  -J  (CM58a)2  +  (CMs5b)2 


with  the  total  trim  angle  calculated  as: 

The  force  trims  are  calculated  in  the  same  manner.  The  trend  of  the  total  trim  angle  from 
the  force  and  moment  calculations  is  shown  in  Fig.  12  versus  Mach  number  and  the 
magnitudes  are  not  the  same  from  the  two  different  methods  at  the  same  Mach  number. 
Usually,  the  trim  moments  are  better  resolved  than  the  normal  force  trims.  In  most  cases, 
the  trim  angles  varied  from  roughly  0.3°  to  a  high  of  0.5°.  The  scatter  in  the  results  is  also 
expected  since  two  bombs  will  not  be  exactly  the  same. 


=  arcsm 
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6.0  COMPARISON  WITH  OTHER  EXPERIMENTAL  RESULTS 

The  MPB-HD  was  extensively  tested  in  the  DREV  indraft  wind  tunnel  (Ref.  14) 
during  the  design  phase  of  the  project  and  as  well  from  full-scale  aircraft  free-flight  trials 
(Ref.  15).  The  MPB-HD  was  also  one  of  the  configurations  that  were  utilized  in  the  Open 
Jet  Facility  (Ref.  2)  that  successfully  showed  that  it  was  possible  to  extract  store 
separation  interference  coefficients  and  free  stream  aerodynamic  coefficients  (static  and 
dynamic),  simultaneously  with  the  MLM  methodology. 

The  reference  center  of  gravity  for  the  moment  coefficient  was  taken  to  be  the 
aeroballistic  range  one,  that  is  3.94  cal  from  the  nose.  The  tabulated  data  from  the  wind 
tunnel,  Open  Jet  Facility  and  from  the  full-scale  free-flight  trials  can  be  found  in  Table 
VIII.  The  aerodynamic  coefficients  measured  from  the  four  experimental  techniques;  the 
aeroballistic  range  (AB-SF,  AB-MF),  the  DREV  indraft  wind  tunnel  (WT),  the  Open  Jet 
Facility  (OJF)  and  the  full-scale  tests  (FF  -  M2470);  are  compared  in  Fig.  13.  The  full- 
scale  tests  only  provided  the  total  drag  coefficient. 

The  axial  force  coefficient  at  zero  angle  of  attack,  CXo  >  is  compared  in  Fig.  13a 

versus  Mach  number.  The  wind  tunnel  data  and  the  Open  Jet  Facility  results  agree 
extremely  well  with  the  aeroballistic  range  data  from  Mach  0.7  to  0.9.  The  full-scale 
results  is  slightly  above  the  aeroballistic  range  results  by  15%.  It  should  be  noted  that  full- 
scale  results  are  for  total  Cq  and  since  the  projectile  flies  at  a  certain  limit  cycle,  it  would 
be  expected  that  the  results  would  be  higher  than  the  Cxo  results. 

CNa,  the  normal  coefficient  slope  versus  Mach  number  is  displayed  in  Fig.  13b 
for  three  of  the  expenmental  techniques.  The  wind  tunnel  and  Open  Jet  Facility  data  is 
within  the  scatter  of  the  aeroballistic  range  results. 
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The  variation  of  the  pitching  moment  coefficient  slope,  CMa ,  with  Mach  number 

is  shown  in  Fig.  13c.  The  lone  Open  Jet  Facility  data  point  as  well  as  the  wind  tunnel 
results  agrees  very  well  with  the  Mach  number  range  of  the  aeroballistic  range  data 


6.1  Limit  Cycle  Amplitude 

As  seen  from  these  aeroballistic  range  results,  the  MPB-HD  is  dynamically 
unstable  at  low  angles  of  attack.  The  motion  plots  (Appendix  A)  show  that,  in  most  cases, 
the  amplitude  at  150.0  m  downrange  is  of  the  order  of  6.0°  to  8.0°  and  still  increasing  but 
leveling  out  in  some  cases.  In  these  aeroballistic  range  trials,  the  initial  angles  of  attack 
were  -relatively  low. 

The  angular  motions  that  were  obtained  in  the  Open  Jet  Facility  tests  (Ref.  2)  at 
roughly  Mach  =  0.75  had  high  initial  angles  of  attack  that  damped  (Fig.  14).  The  angular 
motion  presented  is  only  in  one  plane  and  if  the  projectile  were  slightly  rolling,  that  would 
show  as  low  angles  of  attack  in  one  plane  of  view,  as  in  Fig.  14a.  Nevertheless,  the  motion 
is  definitely  damping  but,  due  to  the  short  time  frame  of  the  tests,  the  final  amplitude  of 
the  angular  motion  was  not  observed. 

From  both  of  these  tests,  the  final  limit  cycle  amplitude  can  be  approximated  to  be 
between  8.0°  to  12.0°.  The  impact  of  this  limit  cycle  amplitude  on  the  range  and 
dispersion  of  the  MPB-HD  is  difficult  to  assess.  Most  of  the  motion  plots  show  a  planar 
motion  at  random  orientations  from  test  to  test,  and  there  is  no  doubt  that  this  will  have  an 
impact  on  the  dispersion  of  this  bomb. 
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7.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  aerodynamic  characteristic  of  a  the  Modular  Practice  Bomb  -  High  Drag 
(MPB-HD)  or  the  BDU-5002/B  MOD  1  were  determined  from  free-flight  tests  conducted 
in  the  DREV  aeroballistic  range.  Eleven  projectiles  were  successfully  fired  in  the  Mach 
number  range  of  0.6  to  0.8.  The  aerodynamic  coefficients  and  stability  derivatives  (Cxo , 

CNcx,  Cjyja ,  and  Cjg5)  were  well  determined.  The  measured  angular  motion  showed  that 
this  projectile  is  dynamically  unstable  at  low  angles  of  attack.  The  pitch  damping 
coefficient  at  zero  angle  of  attack  (CMqo)  and  the  second  order  expansion  term 


(CM  )  were  well  determined.  The  yaw  axial  force  term  as  well  as  pure  side  moments 
qa2 


and  the  trims  were  also  reduced.  A  dynamic  stability  analysis  was  also  conducted. 


Wind  tunnel,  Open  Jet  Facility  experimental  results  as  well  as  full-scale  tests  were 
compared  with  the  aeroballistic  range  results.  The  static  aerodynamic  coefficients  from 
these  three  experimental  techniques  agreed  very  well  with  the  aeroballistic  range  data. 

From  these  aeroballistic  range  tests  and  the  Open  Jet  Facility  experiments,  the  final 
limit  cycle  amplitude  of  the  MPB-HD  can  be  approximated  to  be  between  8.0°  to  12.0°. 
The  impact  of  this  limit  cycle  amplitude  on  the  range  and  dispersion  of  the  MPB-HD  is 
difficult  to  assess.  Most  of  the  motion  plots  show  a  planar  motion  at  random  orientations 
from  test  to  test,  and  there  is  no  doubt  that  this  will  have  an  impact  on  the  dispersion  of 
this  bomb. 


Due  to  the  complex  shape  of  the  MPB-HD,  it  was  not  possible  to  obtain 
aerodynamic  predictions  from  semi-empirical/analytical  codes  since  it  lies  outside  their 
scope  of  application.  It  was  not  deemed  necessary  to  conduct  CFD  calculations  for  this 
configuration. 
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This  database  of  aerodynamic  coefficients  generated  by  this  experimental  program 
is  considered  a  success.  When  it  is  combined  with  the  wind  tunnel,  Open  Jet  Facility  and 
full-scale  free-flight  trials,  it  can  be  used  with  the  BSSM,  the  SSM  as  well  as  the  DREV 
6DOF  trajectory  program  to  predict  accurate  weapon  performance.  Also,  it  can  be 
exploited  with  confidence  for  the  CF-18  OFP.  The  methodology  should  be  extended  to 
other  weapon  systems. 
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TABLE  I 

Nominal  physical  properties  of  model 


d 

(mm) 

m 

(g) 

lx 

(g-cm2) 

(g-cm2) 

1 

(mm) 

CG  from  nose 

(%  /100) 

50.8 

2730.7 

10286.1 

171702.0 

434.34 

0.46 

TABLE  II 

Physical  Drooerties  of  test  projectiles 


Model 

# 

d 

(mm) 

1 

(mm) 

CG  from 

nose 

(mm) 

CG  from 
nose/1 

(-) 

CG  from 

nose 

(cal) 

m 

(g) 

h 

(g  cm2) 

ly 

(g  cm2) 

A01 

50.84 

432.131 

201.809 

0.46701 

3.97 

2750.9 

10296.27 

176578.37 

A02 

50.69 

432.817 

200.792 

0.46392 

3.96 

2709.6 

10141.91 

170918.21 

A03 

50.77 

433.198 

201.317 

0.46472 

3.96 

2721.2 

10213.50 

171677.06 

ACM 

50.77 

433.122 

199.855 

0.46143 

3.94 

2691.1 

10156.59 

166069.14 

A05 

50.88 

433.376 

200.172 

0.46189 

3.93 

2719.4 

10360.04 

167959.99 

A06 

50.85 

433.503 

202.071 

0.46613 

3.97 

2747.3 

10388.09 

173916.72 

A07 

50.80 

433.274 

202.164 

0.46660 

3.98 

2741.3 

10243.65 

174408.93 

A08 

50.85 

432.588 

200.142 

0.46266 

3.94 

2711.2 

10221.22 

168238.05 

A09 

50.93 

432.207 

200.432 

0.46374 

3.94 

2746.6 

10338.28 

172853.03 

A10 

50.88 

432.512 

199.815 

0.46199 

3.93 

2715.9 

10312.85 

168597.78 

All 

50.98 

432.791 

202.072 

0.46690 

3.96 

2778.7 

10425.39 

178226.67 
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TABLE  Vm 
Experimental  results 


a)  Full  scale  free-flight  (Ref.  15) 


Mach 

CD 

0.30 

3.66 

0.40 

3.69 

0.50 

3.72 

0.54 

3.75 

0.60 

3.79 

0.66 

3.84 

0.70 

3.88 

0.74 

3.92 

0.80 

4.01 

0.95 

4.80 

b-)  DREY  Indraft  Wind  Tunnel  (Ref.  U) 


Mach 

CX0 

^Na 

0.5 

3.437 

16.96 

-54.50 

0.6 

3.393 

16.27 

0.7 

3.363 

15.70 

■alia 

0.8 

3.339 

14.78 

-49.16 

0.9 

4.685 

19.02 

-62.40 

1.5 

5.177 

14.27 

-40.23 

c)  DREY  Open  Jet  Facility  (Ref.  2) 


Mach 

Qxo 

cNa 

C-Ma  J 

0.83 

3.379 

18.12 

-48.09  | 
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FIGURE  1  -  MPB-HD  geometry  for  aeroballistic  range  tests 
(all  dimensions  in  caliber,  1  cal  =  50.8  mm) 


t 


I 


JIM 


FIGURE  3  -  Photograph  of  model  and  sabot  package  for  the  MPB-HD  bomb 
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FIGURE  4  -  DREY  aeroballistic  range 
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FIGURE  5  -  Typical  Schlieren  photograph  for  MPB-HD  -  Shot  A 12  (M  =  0.9) 


a)  Station  S24 


b)  Station  S26 
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ARFDAS  -  Aeroballistic  Range  Facility 
Data  Analysis 


Single  & 
Multiple  Fits 


6DOF  Dynamic  Data 

t,x,  y,z,  <t>,  6,  V 


Linear  Theory 
Analysis 

—  — T  ,  _I 

6DOF  Symmetric 
or 

6DOF  Asymmetric 


Fit  Theoretical 
to  Experimental 


Aerodynamic  Forces  &  Moments  vs. 

Mach  No.  &  Angle  of  Attack  &  Roll  Angle 

c  Xo  ,Cxa2,Cxa4’CNa,CNa3’Cf/a5’CYp  a  ,Cyp  a  3  • 
ff3  fl5  mq  >  ^  mq  2*Cnp  a  tCnpcc3>Cnpcc5"- 
Cnp  >  Cnpi ,  C nr  ,  Cnr  2  .  Ca  a  ,  Cm  ya  3  ,  C„  ya  3 
Cip,CfPa2tCts3*Ciycr3 . 


FIGURE  6  -  DREV  Aeroballistic  Range  Facility  Data  Analysis  System 
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FIGURE  7  -  Axial  force  coefficient 


Fig.  7a)  Cxo  versus  Mach  number 


AOA  (deg) 


Fig.  7b)  Cx  versus  angle  of  attack 


Shot  01 
Shot  02 
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FIGURE  8  -  Normal  force  cot 


FIGURE  9  -  Pitch  moment  co< 
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FIGURE  10  -  Pitch  damping  coefficient 


Fig.  10a)  Cjy[qg  versus  Mach  number 


Fig.  10b)  CMq  0  versus  Mach  number 
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APPENDIX  A 


Angular  Motion  Plots 
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